Quantum Phase Transitions in Anti-ferromagnetic Planar Cubic Lattices 
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Motivated by its relation to an A/'P-hard problem, we analyze the ground state properties of anti- 
ferromagnetic Ising-spin networks embedded on planar cubic lattices, under the action of homoge- 
neous transverse and longitudinal magnetic fields. This model exhibits a quantum phase transition 
at critical values of the magnetic field, which can be identified by the entanglement behavior, as well 
as by a Majorization analysis. The scaling of the entanglement in the critical region is in agreement 
with the area law, indicating that even simple systems can support large amounts of quantum cor- 
relations. We study the scaling behavior of low-lying energy gaps for a restricted set of geometries, 
and find that even in this simplified case, it is impossible to predict the asymptotic behavior, with 
the data allowing equally good fits to exponential and power law decays. We can therefore, draw 
no conclusion as to the algorithmic complexity of a quantum adiabatic ground-state search for the 
system. 

PACS numbers: 03.67.-a, 03.65.Ud, 03.67.Lx, 03.67.Mn, 05.50.+q, 05.50.Fh 



Quantum many bodied systems are of increasing inter- 
est in modern physics, particularly systems that exhibit 
a quantum phase transition (QPT). The nature of the 
quantum correlations between the components of these 
systems has been the subject of several recent studies 

QQBBBflftBQGamp.mQ, with the 

suggestion that systems in the vicinity of the critical 
point are highly entangled. In this work we consider 
the phase-structure of a spin-model for which finding the 
ground state can, for certain parameters, be proven to 
belong to the complexity class ./VP-hard 0. We de- 
tect the presence of a QPT between paramagnetic and 
anti-ferromagnetic phases and find that this transition is 
accompanied by a peak in the entanglement between dif- 
ferent parts of the system with a scaling behavior that 
makes it hard to simulate classically [jllhi. Hill. Addi- 
tionally, for the simplest realization of the spin network, 
we have calculated the minimum energy gap for networks 
of up to N = 24 spins. The data accommodates equally 
good fits to both exponential and power law decays, al- 
lowing no conclusion to be formulated as to the efficiency 
of adiabatic quantum algorithms in solving this classi- 
cally A/P-hard problem. Finally, we study the effects of 
frustration in these spin networks, which adds an extra 
element of complexity, altering the phase structure of the 
system. 

The Hamiltonian under consideration is 
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where the sums over and i run over the edges, and 
vertices respectively, of the particular lattice under study. 
The parameters B and V are respectively the longitudinal 
and transverse magnetic fields, and iV is the number of 



spins (qubits), each associated with an individual vertex 
of the lattice. Spin networks of various planar cubic ge- 
ometries are considered 01 1 an d a QPT between param- 
agnetic (large fields), and anti- ferromagnetic (low fields) 
phases is expected in each case. The phase diagram of 
the system is shown schematically in Figfl] 
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Figure 1: Expected phase diagram in the (B,T) plane for the 
systems under study. 

A ground-state search for this system, over the entire 
set of planar cubic lattices, can be shown to be MV- 
hard for the parameters B = 1,T = 0, due to a direct 
mapping to the problem of identifying the maximum car- 
dinality of a stable set, for planar cubic graphs ^5}. In 
the presence of extremely high fields however, the ground 
state is clearly trivial regardless of the size of the system, 
being simply given by the state in which all spins are 
aligned with the field. This gives rise to the possibil- 
ity of performing the, classically ./VP-hard, ground state 
search quantum mechanically, using a quantum adiabatic 
algorithm. The procedure would be as follows: set the 
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quantum register into the paramagnetic ground state of 
the system, with a very large magnetic field, and adiabat- 
ically change this field to the values B = 1,T = 0. Pro- 
vided the evolution is adiabatic, the system is guaranteed 
to evolve to the ground state of this final Hamiltonian 
with a high probability. Clearly the evolution involves 
crossing the phase boundary, from the paramagnetic to 
the anti-ferromagnetic phase, and the time required to 
perform the procedure is determined by the energy spec- 
trum at the critical point. It is, therefore, the scaling 
behavior of the energy spectrum at the critical point, 
that determines the computational cost of the quantum 
adiabatic ground-state search. This provides a strong 
motivation for studying the nature of the QPT in this 
system. 

The set of all possible planar cubic graphs grows 
rapidly with the number of nodes on the graph, with 
the behavior of different instances varying considerably. 
We therefore find it convenient to begin the paper by 
focusing on a particular, simple, instance of the planar 
cubic graphs, which can be easily scaled. We then pro- 
ceed to consider other particular instances, before finally 
considering the average behavior of the entire set. Our 
approach is numerical, and so we are restricted to small 
instances, N < 24 for even the simplest cases, making 
it difficult to draw any conclusions about the asymptotic 
behavior of the system. 

Possibly the simplest subset of the planar cubic lat- 
tices that can be easily scaled, are those consisting of 
two coupled rings, shown in Fig [21 which we refer to as 
the "ladder on a circle" geometry. As with cubic lat- 
tices in general, these lattices always contain an even 
number of nodes N, however, they can be divided into 
two classes. If there is an even number of nodes on each 
ring, the system contains no frustration, while if the rings 
contain an odd number of nodes, then the lattice is frus- 
trated. Although this frustration does not scale with the 
size of the system, the number of frustrated edges in the 
ground state is two, regardless of the system size, and the 
two systems must have identical properties in the ther- 
modynamic (N — ► oo) limit, for the small system sizes 
considered here, the frustration does have some effect. 




Figure 2: "Ladder on a circle" for N = 16 (left) and N = 18 
(right). For N — 16 the system is not frustrated, while for 
N = 18 it exhibits frustration. 



In the case of non-frustrated lattices the phase struc- 
ture of the system is evident even from the small lattice 
sizes considered in this article, N < 24, with clear indi- 
cators of a QPT. Although the frustrated lattice must, 
in the thermodynamic limit, share the same phase dia- 
gram, the phase structure is more difficult to observe in 
small systems. This point is illustrated by considering 
the low-lying energy gaps of the system. Quantum phase 
transitions are marked by the vanishing of the first energy 
gap (A12), the difference between the ground state and 
the first excited state energies. In the "ladder" geometry 
this gap is finite in the paramagnetic phase, decreasing 
abruptly to zero, in the thermodynamic limit, for the 
anti-ferromagnetic phase. In the non-frustrated system, 
the QPT is also marked by a clear minimum in the sec- 
ond energy gap A13, at the critical point, which can be 
easily observed for finite systems, as is shown in FigOfor 
N = 16. For the frustrated ladder however, the multiple 
degeneracy of the ground state in the anti-ferromagnetic 
phase implies that the minimum will only be manifest in 
some higher energy gap, exactly how high is dependent 
of the size of the system, making it far more difficult to 
observe. 

Because of this compression of the energy levels at the 
critical points, it is the energy spectrum at this point that 
determines the minimum time-scale for adiabatic evolu- 
tion between the two phases. In general the adiabatic 
time-scale is determined by the relation 

\(e n (t)\dH(t)/dt\ ei (t))\ 

A2 i ' W 

where |e n (f)) is the n th energy eigenstate of the Hamil- 
tonian H(t). In the case of the non- frustrated "ladder 
on a circle" , the lowest energy level for which the matrix 
element in the numerator does not vanish, is the third en- 
ergy level. In any case, because the ground state in the 
anti-ferromagnetic phase is doubly degenerate, exciting 
the second energy eigenstate during adiabatic evolution 
from the paramagnetic phase does not compromise the 
evolution. It is therefore the second energy gap, the en- 
ergy difference between the ground state and the third 
energy eigenstate, which determines the minimum time- 
scale for adiabatic evolution in this system. The scaling 
of this energy gap with the size of the lattice, then, deter- 
mines the computational complexity of a quantum adia- 
batic ground state search for these lattices. In the case of 
the frustrated ladders however, the multiple degeneracy 
of the anti-ferromagnetic ground state implies that the 
relevant energy gap is higher, exactly how high depends 
on the size of the system, and so it is far more difficult to 
determine the required adiabatic evolution time for such 
a geometry. 

We have analyzed the scaling behavior of the second 
energy gap for non-frustrated ladder geometries for sys- 
tem sizes up to N = 24. The data allow equally good 
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Figure 3: Phase diagram of the non-frustrated "ladder on a 
circle" , for TV = 16 spins. The energy gap between the ground 
and second excited state reveals the appearance of two well 
separated phases. The critical line numerically fits the law 
r = 1.93949 - 0.36730273 + 0.23718273 2 - 0.106107B 3 . 



fits to either exponential or power law decays of the gap 
with lattice size, allowing us to draw no conclusion as to 
the asymptotic behavior of the scaling. A confident de- 
termination of the scaling law in this case would involve 
the analysis of much bigger systems, which are clearly 
outside of our computational capabilities by the direct 
diagonalization techniques employed in this study. As a 
consequence, no conclusion can be formulated as to the 
efficiency of possible adiabatic quantum algorithms for a 
ground-state search for even this restricted subset of pla- 
nar cubic lattices. The problem of computing the quan- 
tum adiabatic complexity of the the classically A/'T'-hard 
ground state of generic anti- ferro mag netic planar cubic 
systems with B = 1 and T = [id l2flj | is even more 
difficult. In this case the difficulty of determining a gen- 
eral scaling law for the energy gap using only a restricted 
number of small lattices is compounded by the fact that it 
is not a priori evident which energy gap determines the 
minimum evolution time. Different instances of planar 
cubic graphs, of the same size, can have vastly differ- 
ent degrees of frustration, and so different ground state 
degeneracies. The energy gap that determines the mini- 
mum adiabatic evolution time will therefore, in general, 
depend on the particular lattice under consideration. 

The two magnetic phases can also be distinguished by 
the entanglement characteristics of their ground states. 
In the paramagnetic phase, the ground state is not en- 
tangled, whereas in the magnetically ordered phase, the 
ground state is highly entangled. The critical point is, 
therefore, marked by some change in the entanglement 
behavior. Indeed it has been argued that one of the sig- 
natures of a QPT is entanglement on all leng th scales at 
the critical point 0, 0, H H H 0, S, IS ES El H El El , 
in rough analogy to classical correlation functions in a 



thermodynamic phase transition. 

In this discussion we use an entanglement measure de- 
fined as the entropy of the reduced density matrix, which 
gives an indication of the entanglement between different 
blocks of particles. For instance, the single particle re- 
duced density matrix obtained by tracing out all but one 
particle, will have an entropy which gives a measure of 
the level of entanglement between that particle and the 
rest of the system. Similarly, if we trace out half of the 
particles, the entropy of the reduced density matrix gives 
a measure of the total entanglement between the block 
of particles remaining, and those that have been traced 
out. 

The ground state entanglement between a single spin 
and the remaining system, calculated in this way for a 
non-frustrated ladder, rises monotonically as the strength 
of the field is decreased, moving from the paramagnetic to 
the ordered phase, due the Zi symmetry of the system 
in the V — > 0, B — > limit (the ground state for the 
"ladder geometry" is close to a GHZ state of N qubits 
for low values of T and B). This degeneracy leads to some 
numerical noise in the data. As any linear combination of 
degenerate eigenstates is also an eigenstate, this implies 
that some quantities, including entanglement measures, 
which differ between the different ground states, become 
ill defined at this point, and depend on which state is 
chosen. Additionally, although actually non-degenerate 
for finite T, in larger systems the first energy gap may 
be too small for the diagonalization algorithm to resolve, 
leading similar problems for small values of T. This noise 
may be seen in some of the plots presented in this article, 
but we point out that wherever it is observed, the limiting 
behavior of the system is clear. 

In contrast to the single particle entropy, the QPT is 
indicated by a peak in the entropy of the reduced density 
matrix of N/2 connected particles S(pm/2)- In the case of 
the ladder geometries we have traced out the particles on 
one of the two rings, and the entanglement calculated in 
this fashion is shown in FigQ] The QPT is therefore iden- 
tified by a collective measure of entanglement, involving 
quantum correlations between large blocks of qubits. 

The scaling of this entropy peak at the critical point 
with the size of the system, seems to be in agreement 
with an area law which predicts that the entanglement 
between two blocks of spins should scale as the size of the 
boundary between the blocks, measured in terms of the 
number of qubits. Note that, despite the fact that our 
configuration is planar, the particular geometry and spin 
partition we have chosen forces the scaling behavior of the 
entropy to be strongly linear, due to the linear number 
of short-range interactions we "cut" with the bipartition 
between the inner and outer ring [1 011 El El El]. This 
result indeed implies that the ground state of the system 
at the critical point is close to a Valence Bond state, as 
described in [1JJ. The scaling of this entropy peak for 
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Figure 4: Entanglement entropy of the reduced density matrix 
for the particles on the inner (outer) ring, for N — 16. The 
peak indicates the crossing of the critical region in the {B,V) 
plane. The critical line numerically fits the law V — 1.93949 — 
0.367302B + 0.237182B 2 - 0.106107B 3 . 



non-frustrated ladders, taken along the line B = 1, indi- 
cates that the scaling tends to obey S{p N / 2 ) ~ 0.08 N 
as the size of the system increases (see Fig|SJ). This be- 
havior is similar to the linear law already found in [j| 
for the Exact Cover adiabatic quantum algorithm and 
clearly differs from the logarithmic law found in [4J, |£( for 
quantum spin chains. Entanglement scaling according 
to this particular bipartition of the system implies that 
even this simple planar model is able to produce an expo- 
nentially large amount of quantum correlations, as mea- 
sured by the maximum rank of the reduced density ma- 
trices over all possible bipartitions. This number can be 
proven to be a measure of entanglement that controls the 
efficiency of certain protocols in order to classically sim- 
ulate the system in a conventional computer. The results 
presented here imply that these classical simulation pro- 
tocols would become inefficient due to the large amount 
of entanglement present in the system || 0, E3, E3 • 

Techniques from Majorization theory EE EH EH 
E3 have proven to be fruitful in characterizing the struc- 
ture of the ground state across a QPT. We have chosen to 
calculate the 2 N cumulants cj = Y^i=i Pi arising from the 
sorted probabilities pj = |(?|V , s)| 2 -where \i) is a quantum 
state, expressed in the z-basis (basis in which the Hamil- 
tonian in EqQ is written), \ip g ) is the ground state of 
the system, and pj are the probabilities sorted into de- 
creasing order- at each step when varying the Hamilto- 
nian parameters. The Majorization analysis of this set of 
probabilities has been useful in the study of efficiency in 
quantum algorithms ESEHEl] Across the QPT on the 
line B = 1, we observe in FigElthat there exists a Ma- 
jorization arrow in the direction of decreasing transverse 
magnetic field T. The phase transition is marked by a 
change in the cumulants at the critical point, reflecting 
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Figure 5: Maximum entanglement entropy per site between 
two rings of a non-frustrated ladder configuration as a func- 
tion of the number of lattice sites, up to N — 24, along the 



line B = 1. For large N the quantity 



*(PiV/2) 



tends to a 



constant, which reveals an asymptotic linear scaling for the 
entropy peak. 
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Figure 6: Majorization arrow for B = 1 when decreasing 
r, for N — 16. There is a change in the behavior of the 
cumulants at the critical point r c ~ 1.8. For simplicity only 
the first five cumulants are plotted. 



the transition from a paramagnetic to a magnetically or- 
dered phase. Step-by-step Majorization implies that the 
ground state is becoming increasingly ordered after each 
step in a very strong sense, at the precise level of each 
one of the probabilities, not just at the global level of 
expectation values such as the mean magnetization. An 
alternative Majorization analysis arises from the proba- 
bility distribution obtained from the eigenvalues of the 
reduced density matrix of N/2 particles, pn/2 EH- In 
the case of the "ladder" geometry we have performed this 
analysis tracing out one of the two rings. The cumulants 
calculated according to this procedure also change their 
behavior at the critical point, this time decreasing as the 
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Figure 7: Majorization analysis of the eigenvalues of the den- 
sity matrix p N / 2 when tracing out one of the rings, when 
decreasing F, for N = 16 and B — 1. Fhere is a change in 
the behavior of the cumulants at the critical point r c ~ 1.8. 
For simplicity only the first five cumulants are plotted. 



m 

Z3 

E 

Z3 
O 



0.8 
0.6 
0.4 
0.2 




□ 



X 

+ +,x 



0.5 1 1.5 2 2.5 3 3.5 4 

r 

Figure 9: Majorization cumulants for the pentagonal lattice 
of N = 20 spins, for B — 1 when decreasing T. A change in 
their behavior is observed at the critical point F c ~ 0.7. For 
simplicity only the first five cumulants are plotted. 




Figure 8: Pentagonal (left) and tetragonal (right) geometries 
for cubic planar lattices of N = 20 and N = 16 spins respec- 
tively. 



field is decreased across the QPT. This can be observed 
in Fig[7| 

Other planar cubic lattices of different geometries show 
similar behavior, with respect to these cumulants, in the 
critical region. For example, we have analyzed the two 
geometries shown in FiglSl (which we call "pentagonal" 
and "tetragonal" geometries) for N = 20 and N = 16 
qubits respectively. These geometries are easily scal- 
able, and for small instances, highly frustrated, although 
the proportion of frustrated edges asymptotes to zero in 
the thermodynamic limit. These lattices also exhibit 
a QPT between paramagnetic and anti-ferromagnetic 
phases, corresponding to a compression of the high en- 
ergy levels of the system at the points at which corre- 
lations, as measured by the entanglement entropy, are 
maximum 10]. 

The Majorization analysis for these lattices reveals 
similar results to the ones presented for the "ladder" ge- 
ometry: for B — 1 and when decreasing T, the cumulants 
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Figure 10: Majorization cumulants for the tetragonal lattice 
of iV = 16 spins, for B = 1 when decreasing V. A change in 
their behavior is observed at the critical point F c ~ 1.6. For 
simplicity only the first five cumulants are plotted. 



obtained from the probabilities for the ground state, in 
the z-basis, increase at the critical point, as is observed 
in Fig[5]and Fig^] In contrast to this, and in agreement 
with the results obtained for the ladder geometries, the 
Majorization behavior of the eigenvalues obtained from 
the reduced density operator pjv/2, decrease at the criti- 
cal point, as is seen in Figlllland Figll2l 

To further illustrate our claims of Majorization signa- 
ture for a QPT, we present results obtained using val- 
ues calculated for all 1249 possible planar cubic graphs 
with N = 18, and averaged. In Fig^|we show the low- 
est two energy gaps, the entropy of the reduced density 
matrix of a single particle (chosen at random), the en- 
tropy of the reduced density matrix of a block of N/2 



6 




1 

0.8 
0.6 
0.4 
0.2 




00 000»o 



■ n + 

□ S + + + 



"12 
A 13 

S[p(2)] 
S[p(N/2)]/N 
C1 Z 
c1„ 



0.5 1 1.5 2 2.5 3 3.5 4 



Figure 11: Majorization analysis of the eigenvalues of the 
density matrix Pm/2 when tracing out half of the system, when 
decreasing F, for N = 20 and B = 1 on the pentagonal lattice. 
Fhere is again a change in the behavior of the cumulants at 
the critical point r c ~ 0.7. For simplicity only the first five 
cumulants are plotted. 
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Figure 12: Majorization analysis of the eigenvalues of the den- 
sity matrix pjv/2 when tracing out half of the system, when 
decreasing F, for N — 16 and B — 1 on the tetragonal lattice. 
There is a change in the behavior of the cumulants at the 
critical point F c ~ 1.6. For simplicity only the first five cu- 
mulants are plotted. Unavoidable numerical noise is present 
for the first cumulant close to F — 1. 

connected spins (chosen at random), and the lowest cu- 
mulant as calculated using the Majorization procedures 
already discussed. The data were taken along the line 
B = 0.5, along which we expect a phase transition in all 
of the geometries. Our results show that all these quan- 
tities change abruptly at the same value of T, strongly 
suggesting the presence of a QPT. 

To summarize, we have studied the quantum phase 
transition in anti-ferromagnetic planar cubic lattices un- 



Figure 13: Here we plot various parameters calculated for all 
1249 possible planar cubic lattices with N = 18, and aver- 
aged. The data were taken along the line B — 0.5, and each 
parameter exhibits a signature of a quantum phase transition. 
The parameters plotted include the first and second energy 
gaps A12, A13, each of which displays a clear minimum. The 
entropy of the reduced density matrix of a single spin (chosen 
at random) S[p(l)] and a block of N/2 connected spins (again 
chosen at random) divided by N S[p(N/2)]/N . Each of these 
quantities shows a peak at the critical region, between V = 1 
and r = 1.4. Finally we have plotted the lowest cumulants of 
a Majorization analysis based on the ordering of the ground 
state in the z basis, cl z , and based on an ordering of the 
eigenvalues of the N/2 particle reduced density matrix cl p , 
The first of these shows a sudden increase, while the second 
exhibits a minimum on the vicinity of the QPT. 



der the action of homogeneous magnetic fields. We have 
shown that techniques from quantum information science 
|30| , such as the behavior of entanglement entropy or Ma- 
jorization theory, can be directly applied to the study of 
quantum many-body systems, leading to new points of 
view in the study of quantum critical phenomena. Lat- 
tices of different geometries display phase diagrams with 
characteristic and differentiated anti-ferromagnetic and 
paramagnetic regions. Scaling laws for entanglement at 
the critical point seem to be in agreement with the area 
law, in a way that even the simplest planar models are 
able to provide exponentially large quantum correlations, 
as measured by the maximum rank of the reduced density 
matrices over all possible bipartitions. The understand- 
ing of the behavior of these systems is well characterized 
by the use of Majorization theory and the analysis of en- 
tanglement entropy, both of which reveal details about 
the transitions present in the models which are some- 
times hard to be obtained by the study of the energy 
spectrum. 
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